Charge instabilities and electron-phonon interaction in the Hubbard-Holstein model 
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We consider the Hubbard-Holstein model in the adiabatic limit to investigate the effects of 
electron-electron interactions on the electron-phonon coupling. To this aim we compute at any 
momentum and filling the static charge susceptibility of the Hubbard model within the Gutzwiller 
approximation and we find that electron-electron correlations effectively screen the electron coupling 
to the lattice. This screening is more effective at large momenta and, as a consequence, the charge- 
density wave phase due to the usual Peierls instability of the Fermi surface momenta is replaced by 
a phase-separation instability when the correlations are sizable. 

PACS numbers: 71.10.Fd,71.45.Lr, 74.72.-h 



I. INTRODUCTION 

In the last years the issue of electron-phonon (e-ph) 
coupling in the presence of strong electron-electron (e- 
e) correlations has been raised in a variety of contexts. 
For instance, in the high-temperature superconducting 
cuprates recent photoemission experiment*^ 3 - indicate 
a sizable coupling of electrons with collective modes, pos- 
sibly of phononic nature, while the softening of a phonon 
peak in inelastic neutron scattering experiments^^ as 
well as features in tunnelling spectral suggest that elec- 
trons are substantially coupled to the lattice in these 
materials. At the same time optical and transport ex- 
periments do not display a strong e-ph coupling ex- 
cept at very small doping, where polaronic features have 
been observed An intriguing interplay between lat- 
tice and electrons has also been invoked to explain trans- 
port in manganites 1 ^, in single-molecule junctions^ and 
in fullereneSfi 2 - where a correlation-enhanced supercon- 
ductivity has also been proposed^ These examples show 
that the issue of e-ph coupling in the presence of strong 
e-e correlation is generally relevant and it translates in 
several related issues. First of all, the fact that vari- 
ous physical quantities appear to be differently affected 
by phonons indicates that the energy and momentum 
structure of the e-ph coupling is important. In turn 
this emphasizes the role of e-e interactions as an effec- 
tive mechanism to induce strong energy and momentum 
dependencies in the e-ph coupling. Secondly, phonons 
may be responsible for charge instabilities. One possi- 
bility is that they mediate interactions between electrons 
on the Fermi surface giving rise to charge-density waves 
or Peierls distortions. It has also been propose d 14 ' 15 that 
a phonon-induced attraction gives rise to an electronic 
phase separation (although in real systems this is ulti- 
mately prevented by the long-ranged Coulombic forces 
with the formation of nano- or mesoscopic domains, the 
so-called frustrated phase separationi 16 ' 17 i 18 : 19 i 20 i 21 i 22 ) 

Due to the above, phonons coupled to strongly cor- 
related electrons have already been investigated by 



means of numerical techniques like quantum Monte 
Carlo^i^^^^ exact diagonalization^^^ Dy- 
namical Mean Field Theory (DMFT), 33,34,35,36,37,38,39 
and (semi)analytical approaches like slave bosons (SB) 
and large-N expansions. 14 i 40 i 41 i 43 i 44 i 45 i 46 Despite this va- 
riety of approaches, a systematic and thorough investi- 
gation within the same technical framework is not yet 
available either due to the demanding character of the 
numerical approaches or to the limited parameter ranges 
investigated so far. Therefore in this paper we study 
the renormalization of the electron-lattice coupling in the 
presence of strong e-e correlations systematically consid- 
ering the momentum, doping, and interaction-strength 
dependencies. In particular we want to elucidate how 
charge density wave (CDW) or phase separation (PS) in- 
stabilities are modified in the presence of e-e interactions. 
To this aim we need a technique which is not numeri- 
cally very demanding, but still provides a quantitatively 
acceptable treatment of the strongly correlated regime. 
In this regard we find the Gutzwiller approach and the 
related Gutzwiller Approximation 47 (GA) a good com- 
promise allowing extensive and systematic exploration of 
various parameter ranges while keeping a reliable treat- 
ment of the low-energy physics. It has recently been 
shown that the Gutzwiller variational approach provides 
remarkably accurate positions of complex magnetic phase 
boundaries in infinite dimensions. 4 - This indicates that 
the Gutzwiller energy and its derivatives are quite accu- 
rate. In this work we extend these results to the charge 
channel also in infinite dimensions, where the GA to the 
Gutzwiller variational problem is exacti 49 ' 50 In addition, 
in order to make contact with layered systems, we study 
the two-dimensional (2d) case where the GA is still ex- 
pected to give an accurate estimate of the energy. This 
also gives us the opportunity to study the interplay be- 
tween nesting in the presence of electron phonon cou- 
pling, which favors Peierls distortions, and strong corre- 
lation which favors phase separation. 

To obtain the phase diagram in the presence of both e- 
e and e-ph interactions, in principle, one should compute 
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the GA energy for every possible charge-ordered state. 
However, it is much more practical to study the static 
response functions of the uniform state to an external 
perturbation and to locate the relevant instabilities. We 
show below that for Holstein phonons in the adiabatic 
limit, the exact charge susceptibility in the presence of e- 
ph and e-e interaction is simply related to the charge sus- 
ceptibility without phonons. Therefore our work reduces 
to compute the latter which is done in the GA. This corre- 
sponds to the static limit of the GA+ (random phase ap- 
proximation), (GA+RPA), which was derived in Ref. l5ll 
and is rooted in Vollhardt's Fermi liquid approach^ As 
a byproduct our work generalizes Vollhardt's computa- 
tion of the zero-momentum and half-filled charge suscep- 
tibility to any momentum and filling. Our approach is 
not as accurate as QMC or DMFT studies as far as the 
electron dynamical excitations are concerned, because it 
inherently deals with the quasiparticlc part of their spec- 
trum. Nonetheless with relatively small numerical efforts 
it allows for a systematic analysis of momentum, doping 
and interaction dependencies of the screening processes 
underlying the quasiparticle charge response and the re- 
lated e-ph coupling. 

It is worth mentioning that the dynamical version of 
our GA+RPA approach has been tested in various situa- 
tions and found to be accurate compared with exact diag- 
onalizatiom 51 ' 53 i 54 ' 55 Computations for realistic models 
have provided a description of different physical quanti- 
ties in accord with experiment , 56 i 57 i 58 The present analy- 
sis of the charge susceptibility of the paramagnetic state 
gives us an opportunity to present the method in a sim- 
pler context with respect to the more complex situations 
considered in the past, thus allowing for the clarification 
of several methodological aspects. On the other hand, 
restricting to the paramagnetic state we ignore antiferro- 
magnetic and related instabilities that arise as the system 
approaches half-filling. 

The scheme of our paper is as follows. We first present 
in Sect. II the derivation of an exact result for the renor- 
malization of the e-ph coupling in the adiabatic limit. 
The GA+RPA approach is presented in Sect. Ill and 
then it is exploited to systematically calculate the mo- 
mentum, doping, and interaction dependencies of the 
charge susceptibility of the Hubbard model. The results 
in d = oo are contained in Sect. IV, while the 2d case 
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where the label "0" indicates that the derivatives must be 
evaluated within the uniform electronic state (in the ab- 
sence of the e-ph coupling); is therefore the ground 
state energy in the absence of the e-ph coupling. The 



is reported in Sect. V. Our conclusions can be found in 
Sect. VI, while the details of our calculations are given 
in the Appendix. 



II. EXACT RELATION BETWEEN 
ELECTRON-PHONON INSTABILITIES AND 
CHARGE SUSCEPTIBILITY 

In this section we show the renormalization of the e-ph 
coupling and how the renormalized coupling is related to 
the electronic susceptibility. We consider a single band 
system with e-e interaction and e-ph coupling on a lattice 

Htot = H e + H ep h + H p h (1) 
Specifically we consider the Holstein interaction 

i i ^ ^ 

where hi is the electronic density operator; n — N p /N is 
the average density on a lattice with N sites and N p par- 
ticles. Pi, Xi and M are the momentum, the displacement 
and the mass of the lattice ions respectively. We treat 
Eq. ^ in the extreme adiabatic limit (M— >oo), where, 
using the Born-Oppenheimer's principle, we can repre- 
sent the ground state of the system as \iptot) — \' l Pe}\Xph}- 
The electronic wave function \ip e ) depends parametrically 
on the ionic displacement Xi . We assume that the ground 
state in the absence of e-ph coupling is uniform [i.e., no 
static CDW state is present]. For a fixed configuration 
of displacements {xi}, the e-ph term of Eq. ^ acts as 
an external field on the electrons, producing a density 
deviation {Sni}, being Srii — {hi) — n. 

The total energy in the adiabatic limit has the form 

E tot = E e [Sn] + E eph [Sn, x] + E ph [x] (3) 

where E e = {^tot\H e \^tot), E eph = (~iptot\H eph \iptot) and 
E p h = (iptot\Hph\iptot) , and Sn, x stand for the sets {cmj 
and {xi} respectively. We move to momentum space and 
perform an expansion of -E^n] up to the second order 
in the density deviation Sn 



(4) 

I 

first order term vanishes identically. For the q ^ terms 
this arises because stability requires that odd powers in 
the <5n q expansion must vanish (otherwise the systems 
would lower its energy by creating a CDW state). In- 
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stead the q = term with dE e /driq = o vanishes because 
we are working at a fixed particle number (<5n q= o = 0). 
Therefore the electronic ground state energy is quadratic 
in the density deviation 



E e [5m] = + 1^2 K^SnqSn- 



(5) 



where we customarily define k^ 1 = ^ gs„2 1 with K q be- 
ing the static charge susceptibility of the electronic sys- 
tem in the absence of the e-ph coupling. 
We now minimize the Holstein energy 



E H = E eph [5n, x] + E ph [x] 



pXiSrii + \kx\ I (6) 



with respect to the ionic displacements Xi at fixed 
5rii finding Xi—~f38rii/ K . Replacing this expression in 
Eq. ([6]) and introducing the adimensional coupling 
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K 



A 



(7) 



[Xo is the density of states (DOS) of the non interacting 
electron system] , we find the following expression for E tot 



Eq. ([12]) has been introduced to separate in the e-ph cou- 
pling the effects of finite q from those of the e-e renormal- 
ization. Specifically, k^/xo contains the effects of finite 
momentum, and is present even for non-interacting elec- 
trons, while the e-e interaction acts on the e-ph coupling 
via a modification of the electronic charge susceptibility, 
given by T q / p. In this exact adiabatic derivation both 
these effects act in a simple multiplicative manner on A. 

A second feature of the result in Eq. §§§ is that the sys- 
tem can become unstable if, upon increasing A, it meets 
the condition A = A c with 



A = A, 



(13) 



for some q c . In this case, if no other (first-order) instabil- 
ities take place before, the system undergoes a transition 
to a charge-ordered state with a typical wavevector q c . 

Now the whole issue to study the effects of e-e inter- 
actions on the e-ph coupling and the related electronic 
charge instability is reduced to the study of the electronic 
charge susceptibility (in the absence of the e-ph interac- 
tion). This is the main goal of the subsequent sections. 



III. FORMALISM 



Em = e£1 + K h ) ~ Sn « Sn -<* (8) 

q 

with 

K CP h = IS: = (Q) 

K " l-XnJx°o 1-A q [) 
where we introduced the renormalized coupling 

A q ^. (10) 

AO 

Eqs. lj8j)- (fT0)) provide the exact second-order expansion 
of the total energy in the adiabatic limit and establish 
a relation between the electronic charge susceptibilities, 
Kq and Kq ph , in the absence and in the presence of the 
e-ph coupling respectively. By construction A q is defined 
in such a way that A q = 1 indicates an instability, in 
analogy with the noninteracting e-ph coupling for which 
A = 1 indicates a q = instability. 

We can write A q = /3r q K q /if, using the charge 
verte x 27 ' 40 ! 43 ' 44 ! 45 ' 46 which acts as a renormalized 
quasiparticle-ph coupling 

r q -> (ii) 

(k is the non interacting susceptibility). Then we find 

Aq ^i-Tq 



In order to compute the static charge susceptibility, 
we evaluate the electronic energy in the presence of an 
external field / 

Hj ^\/-, ; ,.r;.r rT (14) 

ijcr 

where (y i(T {cia) are creation (annihilation) fermionic op- 
erators and fij — f* { . 

Insofar our treatment of the e-e interactions has been 
general. Starting from this section, we will adopt the 
one-band Hubbard model with hopping tij extended to 
whatever neighbors i,j 

H e = Yl tl j c la c ]<y +J2 U ^f ft <l (I 5 ) 

ija i 

fii a =c\ a Ci a is the density operator associated to the op- 
erators cl a , Ci a and U is the on-site Hubbard repulsion. 
In the numerical computations below we consider only 
nearest neighbor hopping tij = —t to be nonzero. 

We apply the Gutzwillcr (GZW) variational method to 
Eq. (TT5]) on d-dimensional hypercubic lattices with lattice 
parameter a=l. We consider the GZW ansatz state \ip) — 
P\Sd), where the GZW projector^ P acts on the Slater 
determinant \Sd). Although we analyze a paramagnetic 
uniform state in order to determine the stability we need 
the energy in the presence of an arbitrary perturbation of 
the charge thus \Sd) shall allow for broken symmetries. 

The GZW variational problem can not be solved ex- 
actly except for particular cases, as for d=oo. Thus one 
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uses the GA. In particular, we use the energy functional 
E e [p, D] obtained by Gebhard^ - which is equivalent to 
the Kotliar-Ruckenstein saddle point energy^ 



with the GA hopping factors Zi„ 



E e [p, D] = ) tijZia-ZjvPjic, + }^UDj 



(16) 



V(i - pa + - A) + y/Di(jm,-<r - A) 

V Piia(l — Piio) 

I 



(17) 



being Pu=^2 a pua- Here p is the single fermion density 

matrix in the uncorrelated state Pj iaa > = {Sd\c\ a Cj a i \Sd) . 

D is the vector of the GA double occupancy parameters 

A = (^|n»T^»il^)- 

To consider arbitrary deviations, the charge and spin 

distribution p and the set of D should be completely un- 
restricted. Since we will consider essentially charge de- 
viations only the part diagonal in the spin indexes con- 
tributes, therefore we will use the notation pij a = pij aa . 

One can show that the expectation value of the diago- 
nal elements pa a , calculated for the \Sd), coincides with 
the value of the density m a , calculated for the GZW state 
Namely, 

n ia = (V'|ct T Cj -|V>) = pua = (Sd\c\ a c la \Sd) (18) 

Eq. (|18p will permit us to express charge density de- 
viations Srii = & n i<y as appearing in Eq. ([3]) via 

We find the saddle point solution minimizing Eq. (|16p 
with respect to p and D. The variation with respect to 
p has to be constrained to the subspace of the Slater 
determinants by imposing the projection condition p=p 2 



6{E e [p,D]-Tr[A(p 2 -p)]}=0 



(19) 



A is the Lagrange parameter matrix. Then it is conve- 
nient to define a GZW Hamiltonian h[p, D] 60 i 61 



hija [p,D] 



dE e 
dp 



(20) 



The variation of Eq. (1191) with respect to p leads to 
h — pA — Ap + A = The Lagrange parameters can be 
eliminated algebraically. 60 Considering also the variation 
with respect to D, we obtain the self-consistent GA equa- 
tions 



[M- = o, 



dE e [p, D] 



0. 



(21) 



(22) 



Eq. (|2"Tj) implies that at the saddle point, h and p can be 
simultaneously diagonalized by a transformation of the 



single fermion orbital basis 



(23) 



leading to a diagonal ho, ho^ — S^e^. Moreover, the 
diagonalized po has an eigenvalue 1 for states below the 
Fermi level (hole states) and for states above the Fermi 
level (particle states). We use a to distinguish quanti- 
ties evaluated at the saddle point. 

In the absence of an external field, we will consider 
the paramagnetic homogeneous state as the saddle point 
solution, i.e. we expand the energy around the param- 
agnetic saddle point, which we describe here. Starting 
from a system with density n = 1 — S (5 is the doping) 
and introducing the notation by Vollhardt, one finds for 
the GA hopping factors^ 



x 



2x 2 



S 2 



1-6 2 



Vl-n + D 

For the GA energy one obtains, 

E e o 
j) 



D. 



Nz^e v + NUD, 
duj wp°(w) 



(24) 
(25) 

(26) 
(27) 



where e , p°, and p denote the energy per site, the den- 
sity of states and the Fermi energy of the non-interacting 
system respectively. 

The minimization of Eq. (|26p yields 



(28) 



x^l-x 2 ) ={l _ 5 2 ) U 



S 2 



8|e°| 



which, by using Eq. (|25p . determines the double- 
occupancy parameter Dq. 

Within the subspace of the Slater determinants, we 
now consider small amplitude deviations of the density 
matrix p due to Hf, given in Eq. IT4]) . This leads to an 
additional contribution Ef[p] to Eq. (JTHJ) 



Ef[p] 



ija 



fijcrPji 



(29) 
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The field / produces small amplitude deviations Sp, SD 
around the unperturbed saddle point density, i.e. Sp = 
p — po and SD = D — Dq: 5p and SD are both linear in 
/. In the presence of the external field /, Eq. (|21[) will 
turn into 



[h + f, P ] = 



(30) 



We will expand Ef A [p, D] = E e [p, D]+Ef[p] around the 
saddle point E e o up to the second order in Sp and J J) 51 i 53 



E? A [p, D] = E e0 + SEP + SEW. 



(31) 



SE^ (SEg) contains first (second) order derivatives of 
the GA energy. We will first consider SE^ in this section 

(2) 

and then SE e ' in the two following subsections. 
The expression for SE^ is 



6EW = Tr[h Sp] + Tr[fSp] 



(32) 



It is convenient to work in the momentum space, where 
^oko-.k'o-' = ^kk"5o-o-' eico- ■ In addition we restrict the ex- 
ternal perturbation to a local field on the charge sector: 
fija — Sijfi with J2i fi — so that Eq. (|29|) becomes 



(33) 



where we introduced the Fourier transform of the density 
deviation <5p q 



1 



fyq = ]y E 5 Pk+iM- (34) 

We will call unoccupied states as particle states (p) and 
use the short hand notation k > \zf for the restriction in 
the momentum; analogously hole (K) state are occupied 
with k < kp, kp being the Fermi momentum. 

The matrix elements of the Sp are not all 
independent^ since p must fulfill the projector condition 
p 2 = p which we can write in terms of Sp 



Sp = p Sp + Sppo + (Sp) 2 



(35) 



Since poko-k'o-' = <5kk'<W'/0k with p^ = 1 for k < kf and 
otherwise po projects onto occupied states. 

Taking matrix elements of Eq. (|35|) one finds for hh 
density deviations (k, k' < kp) 



Sp 



kfjk'cr' 



= E s p 



k(jk"o" 



>Sp 



k"*>k,F ,a" 

and for the pp density deviations (/c, k' > fcp) 

Sp^oKa' = ^ &pVL4k?'<r"$pk**<r"'k*<T t ' (37) 



E 

k" <kp ,<r'' 



Then the hh and pp matrix elements are quadratic in the 
ph and hp Sp matrix elements. Therefore in Eq. (|32p the 
term Tr[ho5p] = Ek e k<5/0kk ^ s nrs ^ or der in the hh and 



pp matrix elements but yields a quadratic contribution in 
the ph and hp matrix elements. The density deviations 
that are off-diagonal in the spin index contribute to the 
magnetic susceptibility^ but not to the charge suscepti- 
bility, therefore in the following they will be neglected. 
One obtains 

Tr[h Sp] = ^2 f kaSpVVa + E £ k'o-<5pk'kv 



k' </cp ,ct 



^ (Cka — ekv)^Pkk'(j^Pk / k<7 



(38) 



k>kp ,k' <.kp;(7 



In the GA the interacting dispersion is related to the 
bare dispersion e£ CT through the relation = z 2 e^ a . 
Eq. (J3H1) shows that the first nonzero contribution be- 
yond the saddle point energy is of second-order in the 
particle-hole density deviations which are our indepen- 
dent variables. 

(2) 

We proceed by considering SEe separately for n = 1 
and general n. 



A. Half-filling case 

Closed formulas can be obtained at half-filling which 
illustrate in a simple manner the physics. This gener- 
alizes the computation done by Vollhardt^ to arbitrary 
momenta. The second order energy contribution for the 
local charge deviations is 



q q 

5nT, u « sd « 5d -* ^ 



being £D q =-^ e-^SDi and 



^q = £ T?(< + +2<_+^_) + 



2(z') 2 
N 



e k+q,o- nko ' 



Q " *2>Z Z rj ^ > q 

L q = 2e z z +D + — — 2^ e k +q,<7"k<7 



(36) n = 2e°z z' c ,+ 



N 



N 



E^ e k+q,CT n ko 



(40) 



where e° = A Eko- e< La n ^- z ' an< ^ z " denote the deriva- 
tives of the hopping factors given in the Appendix. Using 
Eq. ([22]) we can eliminate the SD deviations in Eq. (|39|) 
so that finally the energy functional depends on Sp devi- 
ations alone, i.e. E e JA [p]=E^ A [p,D(p)]. We find 

8D±q = -U„ 1 LqSp±q. 
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Thus the energy E^ [p] is 

kqcr 



with v = 1, ...,d according to the dimension d. Using 
Eqs. (|40]) and (j42]), Aq takes the following form 



where 



q q 



^q = ^q - LlU- 1 



(41) 



(z'l + + 2z'l_+zl_) + 2(z') 2 r ](1 



(42) 



is the GA residual interaction kernel. We have introduced 
the notation Sp^.^, to indicate that only ph elements 
should be taken into account i.e. sums are restricted to 
k > kp and k' < kp. The density deviations can be 
decomposed in ph, hp, pp and hh contributions, ph and 
hp matrix elements contribute quadratically to Eq. (|41|) 
whereas pp and hh are higher order so one should sub- 
stitute 

kcr kcr 

Minimizing Eq. ()41j) with respect to the deviations Sp 
and considering the constraints on the momenta one finds 
the following equation for Spq 



SPq. = -Xq/q ~ Xq^qS p q 



(44) 



Here x q is the static Lindhard function, that is the charge 
susceptibility of the non interacting quasiparticles 



t° = --T 



1 >— ^ ?lk+q,CT — "kcr 



N *— J £k+q : o- — £ko 



(45) 



Notice that ekcr is renormalized by interactions. From 
Eq. (j4"4"|) , we obtain the linear response equatio n 61 i 62 

Spq = -Kqfq (46) 

with the GA+RPA static response function 



Xq 



1 + XqAq 



(47) 



(49) 



From Eq. (|4"9")) we find an analytical expression of the ef- 
fective interac 
U for U < U r 



fective interaction Aq in term of the Coulomb interaction 



_ U(U C + U)(U ~2U C ) 
q AU C (U - U c ) 



(50) 



being U c = 8|e° =1 |. Thus we find that at half-filling 
Aq is independent from the momentum q. In the weak 
coupling limit, we recover the HF-RPA result Aq m U/2. 
Aq is an increasing function of U , diverging at U = U c 



(51) 



Then at the Mott transition not only the charge com- 
pressibility vanishes,— but also the susceptibility K q for 
any momentum [see Eq. Ij47|) ]. 



B. Arbitrary filling case 

The full derivation of SEe is given in the Appendix, 
both in real and momentum space. 

In addition to Spq, we find that it is convenient to 
introduce the quantity 



STi = Uj (5pj i(T + 8p i:ja ) 



and its Fourier transform 

STq = ^(eL + ek+q.J^k+q.crjka (52) 



K q =o is the charge compressibility studied by Vollhardt 
for n = li-S Putting ^4 q =0 and zq = 1, one recovers 
Kq=Xq° for non-interacting electrons. 

In the case of nearest neighbor hopping on a 
d— dimensional cubic lattice we use the relation 



1 \ "* o 

N e k+q, C r"-k 
kcr 



— 22 cos q v = e°riq 



where 



1 d 

Vq = 2 Y C0S q » 
u—1 



(48) 



i5T q corresponds to intersite charge fluctuations, while 
8pq describes the local ones. 

The second order-energy expansion for the charge de- 
viations is given by 



N 



~ ^ Vq5pq8p-q + Zq Z ' d ^ SDqST^q 



^Z (z' + Z ' + _) ^ STqSp^q 

q 

]T LqSpqSD.q + ~ ]T UqSDqSD- 



(53) 
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with the following definitions 



r„ = + 2<_ + z'i_) + ( *' +^- } J2 &K,„"t» 



is the interaction kernel. The elements of IFq are given 
by 



ktr 



4^'+<-) 

TV 



7; £ k+q,CT 



™ko 



krr 



0_« 



2e u z z 



2(^) 2 



^ £ k+q,o- 



«ko 



(54) 



krr 



where z' and z" denote derivatives of the hopping factors 
which are given in the Appendix. 

Using Eq. (f2"2")) one can eliminate the double occupancy 
deviations and arrive at the following functional which 
only depends on the local and intersite charge deviations 



Aq = Vq--± 

Bq = Zo(z' + z' + _) - 2z z' D ^ 



U a 



Aq Bq 

Bq Cq 



where 



Wa 



Aq Bq 
Bq Cq 



Sp-c 

5T- C 



(55) 



(56) 



Since the energy expansion in Eq. (|55|) is a quadratic form 
in 5pq and STq (see also Eq . (f52"j) ) . it is useful to intro- 
duce the following representation for the static Lindhard 
function \q 



N ^ 



c ko- 



"k<7 



£ k+q,(T V c k(T 



(4 C 



fc k+q,(T 
\2 
t k+q,<j) 



^k+q,cr — HjcCT 
^k+q, cr ^kcr 



(57) 



The calculations proceed analogously to the half-filling 
case. We find that the RPA series for the charge exci- 
tations then corresponds to the following Bethe-Salpeter 
equation 



Xq — Xc, 



XqWqXq 



(58) 



For general fillings the response function x q is given by 
a 2 x 2 matrix whose element (xq)n is the charge sus- 
ceptibility Kq. 

For U = oo we can derive an analytical expression for 
Aq valid for any filling n < 1 and any dimension 



Aq = (-e°) 



h — An — Vq 
(1 -n)(2-n) 3 



(59) 



We notice for later use that, since eo < 0, the interaction 
has a minimum at q = 0, a maximum at q = (tt, tt, ...) 
and diverges as n — > 1. Within the GA+RPA the charge 
vertex r q ,2L2242iM«2&<2£ i ntro duced in Section II, is 



given by Zq in the GA. In the following computations for 
r q and g q , we will use = 1. 

In Sections IV and V we will describe the behaviour 
of the charge susceptibility Kq = (xq)n obtained using 
Eq. ([58")) . First we test the performance of the GA+RPA 
approach in d = oo: this is obviously the most suitable 
case for the GA+RPA formalism particularly because the 
GA corresponds to the exact evaluation of the Gutzwiller 
variational wave- function in this limit. We then will move 
to the 2d case, which is physically relevant for quasi 2d- 
materials like cuprates. 

In the following sections we present results for hopping 
restricted to nearest neighbor ty = —t and put t = 1 
which makes the energy and the charge susceptibilities 
dimensionless. Occasionally we will explicitly rescale the 
interaction by the Brinkman-Rice transition U c leaving 
the compressibility units untouched. 



r q = /3 



(Xq) 



(X°q)n 



(60) 



We stress that T q is the renormalized quasiparticle-ph 
coupling. The renormalized coupling for the electrons 
instead is g q = Z q r q , with Zq the quasiparticle weight, 



IV. RESULTS IN INFINITE DIMENSIONS 

We consider the case of a hypercubic lattice in infi- 
nite dimansions with nearest-neighbor hopping, where 



8 



the density of states per spin is given by 



2 1 



7T t 



exp 



2t 2 I 



In this case a momentum dependence in the response is 
still present via the quantity 64 



cos q v 



which enters the interaction kernel Wq and the corre- 
lation functions Xq- For example, the non interacting 
susceptibility reads 




= -4 



duj' 



dw 1 ' 



1 



7T- 



with 



. A q (q/,q;") 

U' - UJ" 

Aq(w>") 
1 



exp 



4z 4 i 2 



(lu'-u") 2 {w' + u"f 



and in the two limiting cases rjq =0 = 1 an( i ^q=Q = ~1 
one can give analytical expressions for the static suscep- 
tibility matrices of the Lindhard function 







\ 2t2 4 



Ei(x) 



where Q is the momentum (±tt, ±7r, ±7r...) and 
denotes the exponential integral^ 

As discussed in Ref. [64], for a generic q the correspond- 
ing 7/ q is trivially zero; only for special q 7/ q takes val- 
ues between -1 and +1. In particular for the hypercubic 
lattice the relevant q's are the ones along the diagonal 
(0,0,0...) — (n, it, 7r...) (and the other equivalent direc- 
tions of the hypercubic lattice, which form a set of mea- 
sure zero). This means that for this infinite dimensional 
lattice, the study of the momentum depencence of the 
quantities is sensitive to effects in the (1,1,1...) direction 
of the Brillouin zone, while it cannot access the other 
directions like, for instance (1,0,0,0...) or (...0,0,0,1). 

In Figs. [1] and [2] we show the dependence on U for 
selected momenta at n — 0.99 (almost half-filling) and 
n = 0.9 respectively. The curves correspond to different 
values of the interaction strength in units of the critical 
U c = 8^/2/7rf, which is the interaction at which in the 
infinite dimensional hypercubic lattice electrons undergo 
the metal-insulator transition at half- filling (n = 1) in 



FIG. 1: d = oo - Charge susceptibility as a function of U/U c 
for n = 0.99. The inset shows an enlargement of the high U 
region. 




FIG. 2: d = oo - Charge susceptibility as a function of U/U c 
for n = 0.9. 



the GA. At small U the susceptibility has a strong en- 
hancement at q c = Q, this is due to the nesting of the 
lattice and leads to the Peierls CDW instability in the 
presence of coupling to the lattice [Eq. (TT5j) ] . 

Starting from the small-C/ side, the charge susceptibil- 
ity is suppressed upon approaching U c and then slightly 
increases again when U is further increased. The behav- 
ior of Kq for U > U c strongly depends on momentum and 
the suppression is most effective for large q and close to 
half-filling. 

Perfect nesting occurs only at half-filling due to the 
matching of the Fermi (hyper)surface when translated 
by Q. In analogy with low dimensional systems one may 
wonder whether away from half-filing an incommensurate 
CDW is favored. 

Figs. [3] and 2] show the momentum dependence of the 
susceptibility for density n = 0.99 and n = 0.9 and vari- 
ous values of U/U c , respectively. 

At small U one finds the nesting-induced enhancement 
at q = Q for both fillings indicating that incommensurate 
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FIG. 3: d — oo - Charge susceptibility as a function of ry q for 
n = 0.99 for various values of the interaction strength U/U c . 
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FIG. 4: d — oo - Charge susceptibility as a function of r? q for 
n = 0.9 for various values of the interaction strength U/U c - 

CDW formation is not favored, i.e. there is no shift of 

Interestingly at large U another instability enters into 
play since in this limit K q acquires a maximum at rjq = 1 . 
One qualitatively recovers a momentum structure simi- 
lar to what is obtained within the large-N expansion of 
the Hubbard model for U = o o 66 i 67 . In this case the 
residual repulsion between quasiparticles is most effec- 
tive at large momenta leading to a suppression of K q for 
r; q — > — 1 . Inclusion of a Holstein coupling would induce a 
PS instability in this limit. The momentum dependence 
of the susceptibility becomes weak for intermediate val- 
ues of U, slightly below U c . All these features are most 
pronounced upon approaching half-filling. 

The fact that for small U/U c the instability momentum 
is pinned at Q is specific to a high-dimensional system, 
where the effects of nesting of the Fermi surface are weak 
and the effects of doping in changing the Fermi surface 
are negligible. We will see in the next section that this 
is not the case in 2d, where upon doping q c moves away 
from the (1,1) direction and shifts along the (1,0) direc- 



tion. 

From Figs. [3] and 2] we conclude that when U is in- 
creased beyond a value of about 0.78C/ C the maximum 
in the charge response moves from large to small mo- 
menta. This signals that the inclusion of (momentum- 
independent) phonons would drive the system towards a 
PS instability at large l/'s, while at small C/'s the systems 
would undergo a transition to a CDW state. The behav- 
ior of the charge response also allows to infer that the e-e 
correlations suppress more severely the e-ph coupling at 
large momentum transfer than at small transferred mo- 
menta. This is the reason why, upon increasing U, the 
system will undergo more easily a low-momentum insta- 
bility (PS), rather than becoming unstable at finite (and 
usually large) momenta. 

V. RESULTS IN TWO DIMENSIONS 

We now move to the 2d case, which is relevant for many 
layered materials. We typically worked with a 100 x 100 
lattice for our computations. 

We start by characterizing the Pcicrls instability in 2d. 
For small U and n — 1 the charge susceptibility has the 
Peierls peak at q = (tt, it), associated with the Fermi 
surface nesting. For the doped system the response ex- 
hibits a peak for a q c close to (ir, 7r), featuring the ten- 
dency to develop an incommensurate CDW in the pres- 
ence of phonons (Fig. [5l full line). The momentum of 
the instability undergoes a shift in the (1, 0) direction of 
the Brillouin zone. For small U the peak is located at 
q c = (a((5)7r, 7r) with 

a(5) = 1 - 0.465 - 1.30S 2 ... ~ 1 - 6/2. 

This depends little on U at weak coupling due to the weak 
q-dependency of Wq close to n = 1. A behaviour com- 
patible to ours has been found also for the 2d-Holstein 
models with QMC and RPA calculations on a 8 x 8 lat- 
tice. 

Along the (1,0) direction « q exhibits another peak at 
q^ (the peak close to q = (0,0), full line in Fig. O. 
This corresponds to the scattering between states at the 
(rounded) corners of the Fermi surface in adjacent Bril- 
louin zones. Upon increasing U (cf. Fig. [5]), the nest- 
ing induced peak structure gets lost. Simultaneously the 
response at large wave-vectors becomes suppressed and 
overcome by that at q = (0,0). This indicates that the 
order of the instabilities is reversed like in the d = oo 
case: For large U the system phase separates before the 
CDW instability arises. This behavior will become more 
clear upon analyzing the charge susceptibility as a func- 
tion of filling and interaction. In Fig. [B] we show the 
charge compressibility in 2d as a function of U/U c with 
U c = I28t/ir 2 . For U = the compressibility is given 
by the noninteracting density of states at the Fermi en- 
ergy. The latter diverges for n — ► 1 due to the Van Hove 
singularity. For U close to U c one recovers a similar be- 
havior as in d — oo. The compressibility vanishes at the 
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(0,0) 



&r.0) 




expansion — . The ground state energy reads 



FIG. 5: 2d-Charge susceptibility for n = 0.9 along the open 
path q = (0,0) - (tt,0) - (tv,tt) for U/U c = 0.5 (solid line), 
U/Uc = 0.9 (dot-dashed line) and U/U c = 2 (dotted line). 
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FIG. 6: 2d-Charge compressibility as a function of U/U c for 
various fillings. 



Mott transition point and has a minimum close to U c for 
n ^ 1. 
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FIG. 7: 2d-Charge compressibility as a function of n for var- 
ious U/Uc- 

In Fig. [7] we present the compressibility as a function 
of n for different values of U/U c . For n — > the com- 
pressibility can be computed exactly using a low density 



iV 



irtn — 2nt 



ln(n) 



where the second term is the leading correction due to 
interactions. Computing the compressibility as = 
d 2 E j dn 2 1 TV one finds that the zero density limit is given 
by the nonintcracting compressibility 

K q=0 («^0) = (27Tt)- 1 . 

In GA+RPA we find instead a small suppression of the 
zero density compressibility with interaction. This is not 
surprising since RPA is expected to break down in the 
low density limit. Still this dependence is quite small 
and we expect that our results are accurate at moderate 
densities. 

For small U the compressibility is an increasing func- 
tion of n and reaches the maximum value for n = 1 as a 
consequence of the van Hove enhancement. For U = U c , 
K q =o goes to zero for n = 1; for larger U > U c , K q =o fiat- 
tens, still exhibiting a smooth maximum for finite doping. 
Therefore the GA compressibility has a jump discontinu- 
ity for n = 1 and U > U c : its left and right limits are 
finite, while it vanishes at n = l42 

The qualitatively different behaviour of the compress- 
ibility for small U and large U is clear: for low fillings 
the system is weakly affected by e-e interactions and its 
compressibility increases with n no matter how large U 
is. Approaching half-filling the correlated nature of the 
system becomes relevant and reduces the compressibil- 
ity of the electron liquid around n = 1. One should also 
keep in mind that close to half-filling AF correlations will 
become relevant. 

In Fig. [5] we show the charge susceptibility /■vq as a 
function of U/U c for n = 1 and n — 0.9 and selected 
momenta. At small momenta the charge susceptibility is 
close to the compressibility in both cases. As the momen- 
tum approaches q = q c the charge susceptibility takes 
its highest values for small U. In particular for n = 1 
(perfect nesting) and U = the charge susceptibility 
diverges indicating that an infinitesimal A renders the 
system unstable. The susceptibility, however, is strongly 
suppressed by U and at half-filling for any momentum 
goes to zero for U = U c . Therefore, as for d — oo, e-e 
interactions renormalize the noninteracting CDW insta- 
bility which thus needs a finite A to occur. 

At small doping K q is finite and shows a shallow min- 
imum close to U = U c . As in d — oo we see from the 
larger K q that for U > U c the PS instability becomes 
dominant. 

The dominance of the PS instability at large U can 
be understood from the strong coupling results close to 
n=l. In this case we expect Eq. (jT7|) . derived for 
n = 1, to be a good approximation. For U 3> U c , K^Aq 
takes large values; then, using Eq. (|4"T|) . we find simply 
that the compressibility saturates as a function of U at 
a value K q « l/j4 q giving a maximum in K q at q = 
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(cf. Eq. (f59|) ). Our results are consistent with Ref. [l5| 
where C. Castellani et al. find PS in a U = 00 slave 
boson (SB) investigation at T = 0. On the other hand, 
SB calculations for 2-d system a 43 i 44 have found PS for 
T t even in the absence of phonons, while for lower T 
an homogeneous state is preferred again. This reentrant 
behavior has not been confirmed by other techniques and 
should be taken with care due to the poor performance 
of mean-field SB techniques at finite T. 

To illustrate the strong coupling behavior, in Fig. [9] we 
show the charge susceptibility for U/U c — 100 along the 
open path q= (0, 0)— (w, 0)— (tt, 7t). If we consider fillings 
quite close to n — 1, the Coulombic repulsion completely 
suppresses the CDW peaks which upon doping become 
visible again for n ~ 0.7. Clearly, more evident CDW 
peaks appear for lower fillings: for n <C 1, the effects 
of e-e interactions are weak even for very large U. For 
n < 0.5, we find that the ground state is a CDW and 
that the susceptibility exhibits quite different features 
with respect to the fillings n > 0.7. 

The most noticeable feature of Fig. [9] is the peak in 
K q centered at q = (0, 0): it gets narrower as the doping 
6 goes to zero. Therefore upon reducing S the charge 
response to a local perturbation spreads more and more 
out in space and the small q peak width is a measure for 
the corresponding inverse screening length. We can give 
an analytical interpretation of this behaviour, using the 
approximate relation K q « 1 /Aq for i5< 1 and adopting 
for Aq its U = 00 form (see Sect. lIIIBl) . We have checked 
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FIG. 8: 2d-Charge susceptibility as a function of U/U c for 
n — 1 (top) and n = 0.9 (bottom). 




FIG. 9: 2d-Charge susceptibility for U/U c = 100 along the 
open path q = (0, 0) — (tt, 0) — (n, n) for n = 0.98 (dashed 
line), n = 0.9 (dotted line), n = 0.7 (solid line) and n = 0.4 
(dotted-dashed line). 
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FIG. 10: Instability surface X(n,U/U c ) of the paramagnet 
towards CDW and PS: the sphered-line marks the transition 
between the CDW state and the state with PS. 



that for low doping the U — 00 form of 1/A q gives a 
quite accurate representation for the whole « q curve. In 
particular, if we consider the low q expansion of l/^4 q , 
we observe that the q = (0,0) peak of the susceptibility 
is well fitted by a Lorentzian peak of half- width q = \/85 
in the x- and y-directions. It is worth noting that the 
typical momentum associated to the peak depends only 
on the doping 8 and not on the energy e . In the low 
q limit a relation q ~ \fS can also be found using the 
U = 00 single SB quasiparticle interaction^ 

All these considerations on the pure electronic response 
lead us to the phase diagram of the e-ph system in 
Fig. QUI The main outcome is that the large e-e inter- 
action changes the nature of the charge instability from 
an incommensurate CDW to a PS. At large U > U c and 
moderate doping this effect already occurs for small val- 
ues of the bare e-ph coupling A. This is the result of a 
compromise between a reduced quasiparticle kinetic en- 



12 



ergy (which renders the system prone to instabilities) and 
the modest screening of the e-ph coupling, when one is 
away from the Mott-insulating phase at half-filling. The 
fact, that the screening is less important at small trans- 
ferred momenta obviously favors the occurrence of PS at 
q c = with respect to the incommensurate CDW. 

In case of the 2-d systems we now consider explicitly 
the renormalized e-ph coupling <? q = z^Tq for the bare 
electrons. In Fig. [IT] we show the behaviour of the quan- 
tity <7 q as a function of the e-e interaction. 

The trend is quite clear: the suppression of the bare 
e-ph coupling is stronger for large q. We compare our 
results at T = with SBi 3 . and QMCSI calculations at 
T = 0.5 (these latters are also performed at a finite Mat- 
subara frequency of the incoming and outgoing fermions 
u> = 7rT). The agreement is generically quite good. How- 
ever, in the finite-T results an upturn of <7 q is also present 
(more pronounced for small q), which was interpreted in 
Ref. |43j as the signature of an incipient PS, which then 
disappears at zero temperature (besides our findings, for 
T = other different treatments find that the ground 
state is homogeneou a 27 i 43 ' 44 ). The nature of this reen- 
trant behavior is still unclear. 

In Fig.rT2]we show the vertex g q for different U along 
a triangular path. This quantity displays minima at 
wavevectors q — q c (and at q^,, see Fig. [S]) thus sup- 
pressing the noninteracting instabilities. These minima 
arise because the bare susceptibility « q is maximal at 
these wavevectors, while the corresponding quantity in 
the presence of U, K q is small due to the suppressed 
scattering at large momenta when the interaction be- 
comes sizable. On the other hand the charge suscepti- 
bility is reduced less at small momenta and this gives 
rise to the pronounced maximum around q=(0,0) in the 
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FIG. 11: 2d- Renormalized e-ph coupling as a function of U/U c 
for n ~ 0.88 at two different momenta. For comparison we 
report the results obtained at finite T (T = 0.5) with SB 4 ^ 
(empty symbols) and QMC? 7 (filled symbols) (these latters 
are for finite Matsubara electron frequency oj = 7r/2). All the 
results are obtained for a 8 x 8 lattice (see Fig. 3 of Ref. HH 
and Fig. 4a of Ref. [13). 



large-?/ case (U > 0.9 in Fig. [T2|) . The shape of the curves 
given in Fig. [T^] is very similar to those obtained within 
a SB calculations at T = 0.00 2. 43 However, this seeming 
agreement has to be taken with a pinch of salt since the 
results shown in Ref. [43| are for the quasiparticle-phonon 
vertex while ours correspond to the vertex for bare elec- 
trons and thus should differ by a factor Zq. 



VI. CONCLUSIONS 

In this work we have investigated the effects of strong 
electronic correlations on the e-ph coupling, in particu- 
lar we considered the case of phonons coupled to the lo- 
cal charge density as described by the Hubbard-Holstein 
model. We first exploited the adiabatic limit of the lat- 
tice degrees of freedom to derive an exact result relating 
the screening of the e-ph coupling to the purely electronic 
static charge susceptibility. This result holds generically 
for any kind of e-e interaction (not only for the Hub- 
bard one) and should also provide valuable information 
in the partially adiabatic case of finite phonon frequency 
(uio « i). It is important to note that the analysis of 
the correlation-driven screening of the e-ph coupling can 
be performed by investigating the purely electronic prob- 
lem. The latter was investigated within the static limit of 
the GA+RPA method. This technique assumes a Fermi- 
liquid ground state and considers the low-energy quasi- 
particle physics. Therefore our low-energy description 
of the electron liquid is appropriate in high dimensions, 
where the Fermi- liquid is a good starting point. Particu- 
larly favorable is the d = oo case, where the GA becomes 
the exact solution of the GZW variational problem. 

The main outcome is that (strong) correlations induce 
rigidity in the charge density fluctuations thereby reduc- 
ing the effective e-ph coupling when this is of the Holstcin 
type. More specifically the analysis of the momentum de- 
pendence shows that the e-ph coupling is more severely 

Q3 gv0 (0,0) (tt,0) {n,n) 



0.3 




(n,7i) (0,0) (tt,0) (n,n) 

q 



FIG. 12: 2d-Renormalized e-ph coupling for n = 0.9 along 
the closed path q = (tt,tv) - (0,0) - (tt,0) for U/U c = 0.46 
(dashed line), U/U c = 0.62 (dotted-dashed line), U/U c = 0.93 
(dotted line) and U/U c = 1.23 (solid line). Results are given 
for a 1000 x 1000 lattice. 
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reduced in processes with large momentum transfer. This 
result, which was already known in large-N approaches to 
the infinite-?/ Hubbard-Holstein models , 14 i 40 i 41 is consid- 
ered here within a systematic variation of the correlation 
strength. In particular, from Figs. QT] and [12] one can 
see that, while at small t/'s the effective e-ph vertices at 
small and large transferred momenta differ at most by 40 
percent, at large f/'s the e-ph coupling at large momenta 
can be five or more times smaller than the couplings at 
low momenta. The fact that the e-ph is screened less for 
small transferred momenta has important consequences 
as far as the charge instabilities of the model are con- 
cerned. Indeed, the e-e interaction not only generically 
reduces the effect of the phonons and enhances the mini- 
mum strength of e-ph coupling to drive the system unsta- 
ble, but also introduces a momentum dependence, which 
changes the nature of the instability upon increasing the 
strength of the correlation: While for small e-e interac- 
tion the leading instability is of the Peierls type, with 
the formation of CDW at momenta |q c | = 2fcp, upon in- 
creasing U, the scattering processes at small transferred 
momentum become comparatively stronger and lead to 
a PS instability at vanishing q c . Our technique allows 
a systematic investigation of how the low-coupling CDW 
instability transforms into the PS instability leading to a 
phase diagram like tha one shown in Fig. \W\ for the 2-d 
system. 

Of course the PS instability is specific to the 

short range nature of the model. When the long 

range Coulomb interaction is included the large-scale 

PS of charged holes is prevented and a frustrated 

PS occurs with the formation of various possible 
texturesi^^o^!^^. 

We also notice that the bare e-ph A needed to drive the 
systems unstable are rather small (of order one or less) 



at large U and moderate doping. A good compromise 
is indeed reached in this region, where the quasiparticles 
have a substantially reduced kinetic energy (the effective 
mass is 3-5 times larger than the bare one), but the sys- 
tem is not too close to the insulating phase, where the 
interaction would screen too severely the e-ph coupling. 
Therefore, in this rather metallic regime the (frustrated) 
PS instability is quite competitive with respect to the po- 
laron formation, which could instead be favored by the 
stronger correlation effects occurring in the antiferromag- 
netic region of the phase diagram^. 

The general interest of the above findings and the en- 
couraging reliability test of the GA technique discussed 
in the present work are a stimulating support for the ex- 
tension of the present work to the dynamical regime. In 
this case future natural extensions will consider the anal- 
ysis of correlation effects on the phonon dynamics and 
their investigation in broken-symmetry states, like the 
stripe phase. 
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APPENDIX 
Real space energy expansion 

Here we give a derivation of the second order term 

SEe in the energy expansion Eq. (|3"Tj) . In real space we 
obtain 



SE^ =J2 f iJ { \z [z {SpuST^ + 6m t STfj) + z' + _{8 Pii 5T? - SrmSTS)] 



+ Z T vo[(z + z + _) 2 5p ii 5p 00 + (z - z + _) 2 5m i 5m j + z Q (z ++ + 2z + _ + z ){Sp u ) 2 

o 

+ z oi z ++ - 2z +- + z )( Sm i) 2 } + \ T ijo[zD{z + z + -)Sp ii SD j + z (z +D + z_ D )5p ii 8D i ] 

+ zoz'nST^SDi + ^oliz'nfSDiSDj + zozUSDi) 2 ] } (61) 

I 

where the fluctuating variables correspond to the which in the charge- and spin sector read as 
charge density p u = p.^ + p^ and the magnetization 
density mi — pu-\ — pu[. Further on we have defined the 
transitive fluctuations 

a 

rp rp + , T C + , T S ^f=E^ C t^)+^W>). 

J-ij - J-tjO + 01 i:j + () 1 , ( a 
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Since we study a paramagnetic system it is convenient 
to define the following abbreviations for the z-factors and 
its derivatives: 

= dz lrT _ , 

Z%<j — o ■ — Z i 

dp~ = Z+ -' dD~ = ZD 

& " & ^ia _ ;/ 3 ^ia " 

d Pu* ~ ++ ' dpiivdpu-a ~ + ~' dp%_ a ~ 

d Zi(j it d Zi(j a d zifj a 

~8D[ = Z °' dpuadDi = Z+D ' dp~~dDi = Z - D 

For the half-filled paramagnetic state we have z = 

and z +D = z_ D . 

Momentum space energy expansion 

We transform Eq. (jrJTj) into momentum space. For the 
paramagnetic system the expansion separates into the 



charge- and spin sector 8E^ — 8E^ + SE^ . 
In the spin sector we find: 

q 

+ ST^SS^) + NJStfS^ ] (62) 
with the following definitions: 

q q 



Aq = [ ( Z - Z +- f e k+q^ n k^ 

The form of SE^ in momentum space is slightly more 
complicate: 



SE * = 2N Z ° {Z ' + Z +-)E^ T -q 

q 

+ Yl^ z ' + z '+-) 2e k+q,a n ^ + {z++ + 2z'+_ + z"_)z e£ ff ?i kCT ]<5p q 5p_ q 

kqcr 

+ ^2 Yy z D( z + z '+-) e i+q,a n ka + z o{ z + D + z'L D )]S pq6 D -q 
kqcr 

+ jv E 2oz D <5L> q 5r c q + E[( z ^) 2e k+ q , CT n k<T + zoz^e^n^SDqSD^ (63) 

q kqcr 

I 

From Eq. (|6"3"|) . we recover Eq. (foil)) by using the defini- tions of L q , J7 q and Vq from Eq. 
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